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We systematically investigate the effect of blockage sites in 
a cellular automaton model for traffic flow. Different schedul- 
ing schemes for the blockage sites are considered. None of 
them returns a linear relationship between the fraction of 
"green" time and the throughput. We use this information 
for a fast implementation of traffic in Dallas. 



I. INTRODUCTION 

In today's crowded world, space and money to build 
transportation systems which can fulfill all demand is 
often not available, or it is not desired to spend it on 
transportation system infrastructure. The result is con- 
gestion: from congested urban centers to congested inner- 
city corridors to congested railways and congested air- 
ports. In consequence, some "forecasting" tool would be 
desirable. Unfortunately, congestion has the side effect 
that causal relations become much more spread both spa- 
tially and temporally M . If a road is crowded, the person 
may attempt a different route or a different mode (spa- 
tial spreading), or she may attempt the trip at a differ- 
ent time (temporal spreading) or even totally drop the 
trip. The result is that planning tools need to consider 
a much wider spatial and temporal context than ever 
before. Conceptually this means that for such problems 
the method needs to be "activity based" , i.e. one needs to 
consider the whole process how people plan transporta- 
tion in a daily or better weekly context (see, e.g., ||). 

Another effect of being in the congested regime is that 
one needs to worry a lot more than before about having 
a dynamically correct representation of the transporta- 
tion system: For example, a peak-period spreading of 
traffic will not show up if one only models a 24-hour av- 
erage situation (as many traditional tools do). Thus, we 
suddenly are faced with a problem where we need to in- 
troduce more dynamical correctness into the modeling 
while at the same time considering much wider temporal 
and spatial scales than before. 

It is fairly obvious that, when faced with a dynami- 
cal problem, a "microscopic" approach, i.e. starting with 
a description of the smallest particles, is in terms of 
methodology the cleanest one. In transportation sci- 
ence, this currently means to consider individual travelers 
rather than, say, aggregated link flows. For example, it is 
difficult to include individual route choice behavior into 
a non-microscopic simulation. There is also some agree- 



ment that the currently most straightforward method to 
deal with microscopic approaches in complicated real- 
world contexts is computer simulation, as opposed to an- 
alytical techniques. Now, when faced with a compute- 
intense problem, such as systematic scenario evaluations 
(see, e.g., or the simulation of the whole national 

transportation system ||, a very detailed and realistic 
microsimulation (see, e.g., ^,^|) may be computationally 
too slow, or too data-intensive to run. 

Alternatives here are simplified models which still cap- 
ture the essentials of the dynamics at the transition to the 
congested regime. Since traffic in general is dominated by 
the bottlenecks in the system, these simulations concen- 
trate on exactly these bottlenecks. The most important 
bottlenecks in urban systems are traffic lights. The nat- 
ural outcome of this way of thinking are queuing-type 
models ||,[l0| • For vehicles that enter the link, one calcu- 
lates when they could arrive at the end of the link. When 
that time is reached in the simulation, they are added to 
a queue at the end of the link. They leave the queue 
once they have advanced to its beginning. The queue 
may have a limited service rate, which models capacity 
restrictions. 

This paper approaches this problem from a slightly dif- 
ferent angle. We use a very simple single-lane microsim- 
ulation to capture at least some of the dynamics that is 
going on on the link itself. This paper will provide a sys- 
tematic approach to such a model. Sec. 2 will describe 
our model, the way capacity restrictions are modeled, 
what their behavior is, and what that means for the re- 
lation between the simulation and reality. In fact, capac- 
ity restrictions are simply modeled by "impurity sites" 
or temporary "blockages" (e.g. Jll|). Sec. 3 discusses an 
implementation and some results for a Dallas study. This 
is followed by a short discussion, highlighting the differ- 
ences between our approach and other "queuing-type" 
approaches (Sec. 4), and a summary. 



II. A SIMPLIFIED APPROACH 

We present a simple simulation model of city traf- 
fic, using a combination of stochastic cellular automata 
(CA) and probabilistic transitions between streets. To 
represent the city network, we use the usual definition 
(e.g. H^l) for links and nodes: a link is a directed street 
segment, such as a bi-directional road divided into two 
links, whereas a node is an intersection; a link can also be 
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defined by an input node and an output node. Vehicles 
are moved on a simple single-lane CA link, and are trans- 
ferred from link to link following a simple probabilistic 
law based on the link's capacity. 



1. Random traffic light 











A. Links 

Links have different characteristics including length, 
speed-limit, number of lanes, maximum capacity, etc. 
The length is necessary to adjust the number of sites 
needed for the discrete approach of the CA. We use 
the standard reference of 7.5 meters for the length of 
one site [Q| 13 1. Each site can be empty, or occupied 
by a vehicle with an integer velocity v 6 {0...v ma x}- 
Vmax — 5 gives good agreement with physical experi- 
ments. 

Since each link is considered as a one-lane segment, 
vehicles are moved using the Nagel/Schreckenberg CA 
rule [l3| . Summarizing the one-lane CA model, the vari- 
able gap gives the number of unoccupied sites in front of 
a vehicle. p no i se is the probability to randomly be slower 
than you could, and rand is a random number between 
and 1. One iteration consists of the following three 
sequential steps which are applied in parallel to all cars: 

1. Acceleration of free vehicles: IF (v < v max ) THEN 
v — v + 1 

2. Slowing down due to other cars: IF (v > gap) 
THEN v=gap 

3. Stochastic driver behavior: IF (v > 0) AND ( rand 
< p nolse ) THEN v = v - 1 

For each link, we introduce an intrinsic probabilistic 
transition, which is a function of the capacity (maximum 
throughput). The one-lane model is faster and easier to 
implement compared to the multi-lane CA. 
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FIG. 1. Schema of the experiment 
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FIG. 2. Flow versus Time for a transition probability of 0.5 



B. Probabilistic transitions 



We introduce various probabilistic models to differen- 
tiate the existing links within a city, from high capacity 
segments such as freeways to low capacity segments such 
as arterials. If we consider only one-lane links, the proba- 
bilistic transition is introduced to control the output flow 
of a link. A high capacity link will produce a high output 
flow, while a low capacity link will produce a low output 
flow. 
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FIG. 5. Space time diagram for p=0.2 



Let us consider the experiment in fig|l|], consisting of 
two consecutive links separated by a probabilistic transi- 
tion ptrans- The first site of link 1 operates as a generator 
of vehicles, where one vehicle is introduced per n itera- 
tion^). The flow measured at the end of the second link 
versus the number of iterations is shown in fig^). The 
probabilistic transition is set to 0.5 in this example. The 
flow measured at iteration t is the number of vehicles that 
left the second link until that moment, divided by t. As 
a result, the unity of the flow is vehicle per iteration. We 
introduce 1 vehicle every 3 iterations at the first site of 
the first link with maximum velocity 5. This is enough to 
assure that the first link will reach around 1200 veh/h for 
a Pnoise of 0.5, which is close to the maximum through- 
put of such a link in the CA implementation O] . If the 
first site is not empty at the introducing time step, we do 
not add vehicles. The vehicle's velocities are updated by 
the one-lane CA model before reaching the intersection. 
If the vehicle is allowed to go through the intersection by 
the CA forward rule, we check the probabilistic transi- 
tion. 

If the generated random number is lower than the 
probability ptrans, the vehicle keeps its velocity and 



reaches the second link. In contrast, if the random num- 
ber is greater than ptrans, we place a fictitious car in the 
first site of the second link in order to force the vehicle to 
brake and stop at the intersection. Technically: If a car 
reaches the last five sites of a link, it produces a random 
number. We introduce the simple algorithm: 

1. Transition check: 

IF (rand < Ptrans) THEN normal C A- update 
ELSE gap=distance from the vehicle to the inter- 
section 

This situation is in principle well understood. The "im- 
purity site" will create a reduced flow that can pass that 
site, and since flow needs to be conserved along the link, 
this sets the maximum throughput for the link [|lT| , p^ — p~7| . 
Yet, in the context of the stochastic traffic cellular au- 
tomaton as used here, we are only aware of Ref. Jl8| , 
and the specific mechanism used there is not the one we 
wanted to use. 

Figs |||4|] demonstrate the formation of traffic jams 
spreading to the beginning of the link, caused by braking 
of vehicles. The beginning of the second link can again 
be considered as a generator of vehicles. Nevertheless, 
the input flow and p are not proportional. 

To illustrate this comment, we conduct the same ex- 
periment with probabilistic transitions ranging from to 
1. The average flow obtained for each experiment is pre- 
sented in figP]. For each data point, the flow is averaged 
in the time period (5000,20000). See fig§. The intersec- 
tion does not function as a perfect generator of service 
rate p. If a vehicle leaves the last site of the first link, this 
vehicle is not automatically replaced, due to the stochas- 
tic third step included in the one-lane CA model. This 
plot can be divided into three different parts: 




Transition Probability 
FIG. 6. Flow versus transition probabilities 

(i) A high probabilistic transition {ptrans between 0.8 
and 1.0) gives linear results with input flow. In this sce- 
nario, vehicles do not stop often at the intersection, thus 
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the intersection does not work like a stop and start point. 
See fig§. 

(ii) A low probabilistic transition (between and 0.4) 
gives results that can be explained by a simple hypothe- 
sis. Most of the cars stop at the intersection and form a 
compact traffic jam, as shown in fi.g||. There are no im- 
portant spaces in this queue. Assuming that the second 
last site of the first link is always crowded, how many 
iterations does a vehicle need to go through the intersec- 
tion? If a vehicle is on the last site of link 1, the vehi- 
cle needs l/p no ise iterations on average to advance, and 
then multiplied by 1/ptrans to go through the intersec- 
tion. Viewed from the perspective of the next following 
vehicle, that one needs to wait l/p no ise ■ ^/ Ptrans steps 
until the vehicle ahead is gone, and then another l/p no ise 
steps to move itself to arrive at the last site. As a result 
the average number of iterations for a vehicle to advance 
from the second last site of link 1, to the intersection is 
l/Pnoise + l/(PnoisePtrans)- This could in theory be con- 
tinued, but it would not necessarily get better because 
one would need to include the influence of "holes" in the 
queue; or, more technically: The approximation is only 
valid for ptrans — ► 0, and second order corrections are 
thus negligible. In any case, the corresponding flow is 

F ~ PnoisePtrans / (1 ~t~ Ptrans) ■ 

The function, F, shown in fig@ fits well to the data 
measured for low values of p, while for p > 0.4 the hy- 
pothesis is no longer valid. 

(iii) Figure || demonstrates what happens for proba- 
bilistic transitions between 0.4 and 0.8 at a microscopic 
level. Within the queue, holes are generated by the inter- 
section and an analytical approach becomes more diffi- 
cult. Periodically, vehicles pass through the intersection 
without braking and stopping, which produces a higher 
flow compared to the linear relationship illustrated in 
ng§- 

Many experiments can be conducted using other prob- 
ability distributions for the intersection. The model pre- 
viously described operates like a random traffic light, 
where the light becomes green with the probability 
Ptrans, which is also the fraction of the time the light 
is green: f gre en = Ptrans- This model can be con- 
sidered to be one between two extreme distributions, 
where in between the extreme cases one can encounter 
an infinite number of distributions that keep the frac- 
tion of a green light of the total time of a traffic cy- 
cle constant. The first distribution is a classical traf- 
fic light. The green fraction here is straightforward: 
fgreen = T green / (T green +T red ) . We call the second model 
a Dirac traffic light. As we work with discrete systems, 
the objective is to set a green light or a red light on only 
one time unit, equally spaced on a cycle. The green frac- 
tion is fgreen = l/(l+T red ) for T red > 1 (and T green = 1 
by definition) or I — 1/(1 + T green ) for T green > I. All 
three distributions are illustrated in fig . 
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FIG. 7. Different distribution probabilities 

Next, we present the same experiments discussed 
above, for these two distributions. 



2. Normal traffic light 

We repeat the same experiment described in figQ| with 
a normal traffic light at the intersection. 
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FIG. 8. Space-time diagram normal traffic light p=0.5 FIG. 9. Space-time diagram normal traffic light p=0.9 

The dissolution of a queue as the light turns period- 
ically green is shown in ngQ. This phenomenon does 
not provide an easy analytical solution. For each green 
fraction / greerl ranging from to 1, the input flow of the 
second link is measured and is illustrated in ngjnj. This 
relationship is almost linear. For high values of transi- 
tions, vehicles still have to stop occasionally, which de- 
creases the output flow. Figure l9|, when compared to the 
space-time diagram produced by the random traffic light 
figQ, displays a lack of fluidity. 
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FIG. 10. Flow versus Transition probability 



3. Dirac traffic light 

The Dirac traffic light generates the highest flow for a 
given f g ree n in the experiment or figjfj. The space-time 
diagram performed with a probabilistic transition of 0.5 
is given as an example. In this case, the traffic light is 
successively green and red. Figure |ll]] shows less com- 
pact traffic jams at the end of the first link than the other 
space-time diagrams for the same probabilistic transition. 
This is still due to the vehicles that pass through the in- 
tersection at maximum velocity without braking. The 
analytic explanation for this is the fact that the paral- 
lel update tends to generate states where particles are 
followed by holes, sometimes called "particle-hole attrac- 
tion" Hi. 




FIG. 11. Space-time diagram Dirac traffic light p=0.5 
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Dirac traffic light p=0.16 



FIG. 13. Space-time diagram Dirac traffic light p=0.9 

The output flow of link 1 for any value of p is much 
higher than the two flows measured previously for the 
two other probability distributions. There is no linear 
relation at any position on this diagram. The space-time 
diagrams plotted for a p = 0.16 and p — 0.9 exhibit more 
fluidity for the output traffic fig p2|Jl^ ]. 
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Transiton probability 
FIG. 14. Flow versus transition probability 

III. DALLAS 
A. Implementation 

The normal traffic light model is the most linear model 
simulated in this paper. On the other hand, setting a 
traffic light at each intersection would cost computation 
time. The random traffic light presents the advantage to 
be checked only when a vehicle reaches the intersection. 
The vehicle generates a random number which allows it 
to drive trough the intersection or not. 

We apply this model to the Dallas Fort- Worth area. 
The context is the so-called Dallas/Fort Worth case 
study |4|j20|] which has been done as part of the TRAN- 
SIMS (Transportation ANalysis and SIMulation Sys- 
tem) project gj. TRANSIMS uses individual route plans 
for each individual traveler. A route plan consists of a 
starting time, a starting location, a list of links the ve- 
hicle intends to follow, and an ending location. A mi- 
crosimulation in the TRANSIMS project such as the one 
described here is thus faced with the task to move these 
vehicles according to these specifications. 

One immediately observes that one somehow has to 
correct for the fact that we are only using single-lane 
roads, i.e. our links will usually not be able to carry the 
prescribed number of vehicles. We solve that problem by 
using a sub-sample of the plans. The size of that sub- 
sample is obtained as follows: 

• Pnoise = 0.5 results in a maximum throughput of a 
link of approximately 1200 veh/h (using ptrans = 
1). 

• We search for the link with the highest capacity in 
the area we want to simulate. In our case, this was 
a four lane freeway with a capacity of 7800 veh/h. 

• We thus need to sub-sample the population by a 
factor of 1200/7800 w 0.154, i.e. a route plan from 



the full plan-set is going to be used with a proba- 
bility of 0.154. 

• Links which have a lower capacity than 7800 veh/h 
take this fact into account by using a value of ptrans 
according to Fig. [| i.e. if the value of the road is 
C, then the value C ■ 0.154 is used on the y-axis to 
find the correct value of ptrans on the x-axis. 

A more precise calibration is more complicated than this 
because it also depends on the interplay between route 
planning and route execution. This is clearly out of the 
scope of this paper; further publications on the subject 
are in preparation. 

B. Simulation results 

In this section, we want to give some examples how 
this simulation is going to be used. These examples will 
be given in the context of the TRANSIMS Dallas/Fort 
Worth case study. That case study used as input a street 
network of the Dallas/Fort Worth area, containing 24662 
links and 9864 nodes, and information of all trips in 
this area during a 24 hour period (approx. 10 million 
trips). The study focused on a busy 5 miles times 5 miles 
area north of downtown Dallas, and on the time between 
5am and 10am. This still involved 300000 trips. As 
mentioned above, micro-simulations in the TRANSIMS 
project are route-plan driven. Thus, for each of these 
300 000 trips, route plans were calculated. The fact that 
drivers adjust to congestion caused by other drivers was 
taken into account by iteration several times between the 
route planning and the micro-simulation. For further in- 
formation, see Refs. ]2(i|-|22]] . 

One important specification missing in the above de- 
scription of the micro-simulation is how vehicles enter 
and leave the simulation. TRANSIMS specifies parking 
locations along links, which represent all parking oppor- 
tunities that can be reached from this link. In order to 
prevent that the traffic that leaves parking unrealistically 
disturbs the traffic flow, vehicles from the parking loca- 
tions are only inserted if v max sites backwards from the 
parking location are empty. If the space is not free, the 
car is placed in a queue, waiting to enter the simulation 
in one of the following iterations. 
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FIG. 15. Snapshot of the case study at 7:00 

A snapshot of such a simulation with the model de- 
scribed in this paper can be found in Fig[fil|. The denser 
square area in the center represents the study area, where 
all streets including local streets were represented in the 
data base. For this example, also the streets outside that 
area were simulated. Dots denote individual vehicles. In 
this plot, most of the traffic is on the freeways, as is real- 
istic. Also, one notes that for lower capacity road, traffic 
is mostly queued up towards the end, as one would expect 
from the dynamics of the model. Yet, this is really not 
too unrealistic since also in reality traffic through minor 
roads tends to queue up at the ends. 
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FIG. 16. Space-time plot of a particular link (Beltline Rd., 
an east-west arterial in the northern part of the area) from 
7:00am to 7:05am (left), 7:30am to 7:35am (middle), and 
8:00am to 8:05am (right.) 

The space-time diagram of five consecutive links is 
shown in figjlq]. These links are a part of an east- west 
arterial located in the north of the study area. The figure 
shows nicely how queues built up at the end of links due 
to the capacity restrictions. 
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FIG. 17. Travel times versus departure time 
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Many statistics can be extracted from the simulation. 
As a further example, we present the travel time versus 
departure time for each vehicle (Fig. |l^). This figure 
shows nicely that even such a simple simulation as the 
one described in this paper can, given a realistic trip 
demand input, display the higher travel times during the 
rush hour. 
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C. Computational performance 
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FIG. 18. Real time ratio 

We present a performance diagram in fig|Tg|] where we 
introduce the RTR versus the simulation time. The RTR 
is the ratio of the real time on the simulation time. This 
example of simulation was executed on a SUN UltraSparc 
CPU with 250 MHz where approximately 46000 plans 
were simulated in the whole Dallas Fort- Worth area. The 
diagram ngp8| shows a ratio of 23 in the middle of the 
rush period, but in average the ratio is around 28. This 
clearly shows that simulations like the one described here 
have enough computational speed for thorough investiga- 
tions of traffic problems. 

IV. DISCUSSION 

Transportation models using simplified link dynamics 
ultimately fail to generate some aspects of a complicated 
reality, such as a turn pocket having a queue spill-back 
into the lanes going straight. Yet, as pointed out in the 
introduction, using a highly realistic model sometimes 
is not an option, for example because of computational 
constraints or data collection constraints. In such cases, 
knowing the different limitations of the simplified models 
becomes crucial. 

For example, the different simplified models handle 
queue discharge dynamics in different ways. In our 
model, when a vehicle is moved to the next link, this 
leaves a "hole" on the link where the vehicle comes from. 



In the next time step, that hole may or may not be filled 
by an advancing vehicle, according to the stochastic driv- 
ing rules. In a congested situation, this hole slowly trav- 
els backwards, until it eventually reaches the other end 
of the link, allowing a new vehicle to enter the link. In 
Ref. pQl , holes are transmitted instantaneously to the 
other end of the link. The method of Ref. || assumes 
infinite queuing capacity on each link. It is clear that all 
three methods will generate different dynamics. 

As another example, in fairly realistic models, sources 
and sinks for traffic are better located in the middle of 
links instead of at nodes. Vehicles attempt to squeeze in 
between other vehicles at that location on the link. If 
the link gets congested, the additional vehicles will have 
trouble finding additional space to squeeze in. It is clear 
that models who totally give up a representation of traffic 
dynamics on the link such as Jl0|,|| will lead to different 
behavior for traffic sources and sinks. 

Certainly, the simple queuing models could compen- 
sate for that. Yet, that usually comes at the price of being 
tedious. Often, it will be more straightforward to move 
directly to a higher fidelity (but usually computationally 
slower and more data intensive) micro-simulation. We 
believe that, at the current stage, it is more important 
to really understand the dynamical differences between 
different models and to compare their behavior in real- 
world applications, than to attempt to improve simplified 
models in non-intuitive ways. 

Last but not least, the model presented in this paper 
actively moves vehicles along links with roughly realistic 
dynamics. This makes graphical output such as in Fig. [l5| 
much more intuitive and appealing. 

V. SUMMARY 

"Blockage" sites, i.e. sites which move particles or ve- 
hicles only a fraction of the time, reduce the maximum 
throughput of a link of cellular automata models for traf- 
fic flow and particle movement studies. We have sys- 
tematically tested the effects of three different blockage 
schemes, where one was the usual random draw, one was 
a regular traffic light with long red and green times, and 
one was what we called a "Dirac" traffic light because 
it had 1-second spikes of one color. In general, there is 
no linear relation between the fraction of green time and 
the throughput. The Dirac traffic light returned the high- 
est throughput; the explanation for this is the "particle- 
hole" attraction that can be found in the type of cellu- 
lar automaton that was used. Since none of the timing 
schemes returns a totally linear relation, we used the ran- 
dom scheme in an implementation of traffic in Dallas. We 
showed some exemplary results of this implementation. 
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